arXiv:1503.08507vl [math.DS] 29 Mar 2015 


On a planar circular restricted charged three-body 
problem 

Abimael Bengochea • Claudio Vidal 


Abstract We introduce a circular restricted charged 
three-body problem on the plane. In this model, the 
gravitational and Coulomb forces, due to the primary 
bodies, act on a test particle; the net force exerted by 
some primary body on the test particle can be attrac¬ 
tive, repulsive or null. The restricted problem is ob¬ 
tained by the general planar charged three-body prob¬ 
lem considering one mass of the three bodies going to 
zero. We obtain necessary restrictions for the param¬ 
eters that appear in the problem, in order to be well 
defined. Taking into account such restriction, we study 
the existence and linear stability of the triangular equi¬ 
librium solutions, as well as its location in the config¬ 
uration space. We also obtain necessary and sufficient 
conditions for the existence of the collinear equilibrium 
solutions. 

Keywords Hamiltonian Vector Fields, charged three- 
body problem, circular restricted charged three-body 
problem. 


1 Introduction 


The force exerted by electric E and magnetic field B on 
a particle of velocity v and charge g is F = dE+vxB), 
and it is deno minated Lorentz’s force ( Goldstein 198d 
Jackson 1975h . The exact description of this problem 


is formulated in the relativistic context. Nevertheless, 
when the speed of the particle is much smaller than 
the one of the light c (around 300, OOOkm/s), the dom¬ 
inant term is the Coulombian since the first relativistic 
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correction is of order (u/c) 2 , as it is shown by the La- 
grangian of Darwin ( Jack son 19751) . In the study of 
charged particles, the gravitational interaction is com¬ 
monly of smaller magnitude than the electromagnetic 
one, whereas in celestial problems it happens the oppo¬ 
site. 

Several works have been written about restricted 
models concerning three point particles (it is assumed 
that one of the bodies does not affect the motion of 
the other two particles) with additional interactions to 
the gravitational one, for instance the Coulomb or pho- 
togravitational cases (the gravitational force is attrac¬ 
tive, the Coulomb can be either attractive or repulsive, 
and the photogravitational only repulsive, in a generic 
w ay - the zero value is allo wed in all of them). 

RadzievskiH ( 195(1 19531 ) introduced the photogravi 


tation al restricted three-body problem, and lDion vsiou a nd Stam oul 
1989 ) a restricted charged three-body problem. ’Scliuermanfr 


198011 studied the stability, and location on the config¬ 


uration space, of the equilibrium solutions when radia¬ 
tion pressure and Poynting-Ro berts forces ar e inclu ded. 

In the photogravitational case, Kunitsvn and Tureshbaev^ 
19831 ) studied the exist ence and stability of the collinearj 


equilibrium solutions, Lukvanov (198 41 the existence of 
the triangular equilibrium solutions, and the collinear 
equilibriu m so lution, in a paramet ric way. In the same 
model, Kunits yn and. Tu reshbaev (ll985aT) considered 
the existe nce o f t he co llin ear equilib rium solutions, 
and Kunitsvn and Tureshbaev ( 1985bh the existence 
and stab ility of the equilibrium triangular solutions. 
Simmons ( 1985t l studied the restricted case with radi¬ 


ation pressure; part of his study is concerned with the 
analysis of the triangular, collinear and spatial equi¬ 
librium solutions. In the photogravitational restricted 
problem, Lukyanov (|1986n studied the stability of the 
collinear and triangular equilibrium solutions, as well 
as its location in the configuration space. 
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Some of the studies present the same potential func¬ 
tion, including our own. However, the parameters, 
as well as its allowed values, can be different. For 
insta nce, in the restricted char ged t hree-body prob¬ 
lem, Dionysiou and Vaiopoulps (119871 ) consider a dif¬ 
ferent parameterization from the one used in this work. 
In fact, the potential function is V = 
where p £ ( 0 , 1 / 2 ], q £ R, and p, denotes the dis¬ 
tance between the body i = 1,2 and the third body. 
On the other han d , in the photogravitation a l case , 
Radzievskii ( 195(ih . Kunitsyn and Tureshbaevl ( 1983 ). 
Lukvanov (1984) and Simmons ( 19851) . consider a pa¬ 
rameterization similar to ours. The potential function 
is V = 71 where p, 71 and 72 are real param¬ 

eters, and pi, i = 1 , 2 , the distance between the bodies, 
as in the charged case. The parameter p belongs to 
(0,1/2], whereas 7 i £ (— 00 ,1], i = 1,2. The qualita¬ 
tive properties of the gravitational and photogravita¬ 
tional forces acting on the test body, by the primaries 
i = 1 , 2 , are defined by the parameter 7 ,: if 7 , < 0 
the gravitational force dominates over the photogravi¬ 
tational one, for 7 * = 0 both forces are equal to each 
other in magnitude, for 0 < 7 * < 1 the gravitational is 
the strongest one, and for 7 * = 1 the photogravitational 
force is zero. The case 7 * > 1 has not physical sense 
because it corresponds to an attractive photogravita¬ 
tional force; usually it is included in the study of the 
mathematical model. The different selection of units 
and generality with which these problems have been 
studied can be appreciated in the number of parame¬ 
ters that appear, that is, two and three respectively. 

This paper is organized as follows. In Section 2 we 
give a brief review of the problem of two charged bod¬ 
ies. Here we emphasize that the necessary condition 
on the parameters of mass to, and charge qi, i = 1 , 2 , 
for the existence of circular orbits of the particles 1 
and 2, is G — k > 0 (G is the constant of uni¬ 

versal gravitation and k is the constant of Coulomb); 
this limits the values of the parameters that appear 
in the restricted problem of three charged bodies on 
the plane, which we enunciate in Section 3. With 
this aim, we make the mass of the third body tends 
to zero, which gives us a well defined restricted circu¬ 
lar problem. The potential function associated to this 
problem is given by V = ^ , with param¬ 

eters p, /?i, /?2 , and the usual distances p\, p 2 - Here 
p £ (0, 1 / 2 ] and fa £ R for i = 1 , 2 , with the restriction 
(/?i — l )(/?2 — 1 ) < 1 , which is the necessary condi¬ 
tion for having circular solutions for the bodies 1 and 
2. Such condition has not been considered in previous 


studies, for instance Dionysiou and Vaioooulos (1987) 


Dionvsiou and Stamoul dl989l).lKunitsvn and Tureshbaevl 


( 19831) . Kunits yn and Tureshbaev (|l985b ). L ukvan ov 


(1984), Lukvanov ( 1986 ). Simmond ( 1985 ). reason why 
our problem acquires a great difference to others al¬ 
ready treated and justifies its study. We remark that, 
if the parameters that appear in the potential are not 
properly considered, the interaction between the pri¬ 
mary bodies will be repulsive and might prevent those 
bodies from having a circular movement. As in the pho¬ 
togravitational model, the parameters i = 1,2 de¬ 
termine the qualitative features of the gravitational and 
Coulomb forces. Actually, we have a similar descrip¬ 
tion for both photogravitational and charged problems, 
with the exception that /3, > 1 is physically possible 
and corresponds to a Coulomb attractive force. Once 
established the restricted charged three-body problem, 
in Section 4 we study the existence of triangular and 
collinear equilibrium solutions. In Section 5 we deal 
with the linear stability of the triangular equilibrium 
solutions and its location in the configuration and pa¬ 
rameters space. We finish the article with the conclu¬ 
sions of this work. 


2 Dynamics of the two-charged problem 


Consider an inertial frame of reference. The Hamilto¬ 
nian associated to bodies 1 and 2 of charges q\ and q 2 , 
and masses m\ and m 2 respectively, with gravitational 
and Coulombian interaction, is 


U(ri,r 2 ,pi,p 2 ) 


IIPill 2 IIP2II 2 kqiq 2 - Gm 1 m 2 
2m 1 2 m 2 ||ri — r 2 || 

(1) 


where k > 0 is the Coulomb constant, G > 0 is the 
universal gravitational constant, and 17 £ R 3 , p, £ R 3 , 
i = 1,2 are the positions and moments of the bodies 
with the same index, respectively. 

Given the Hamiltonian ©, and the initial condi¬ 
tions, the equations of Hamilton f ,. = H Pi , p, = — H ri 
determine the motion of each body. The equations of 
motion correspond to a system of differential equations 
of second order: 


TOlfi 


— (Gmim 2 - kqiq 2 ) 77 —^ 

||n — r 2 || 3 


m 2 f 2 = — (Gmim 2 — kqiq 2 ) 


r2 - ri 

Il r i — r 2 1| 3 ’ 


where the dot denotes derivative with respect to time 
t. In order to classify the solutions, we define C = 
Gmim .2 — kq\q 2 - Three cases are identified: 


• C> 0, 

• <7 = 0, 

• c < 0 . 
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The first one is equivalent to the Kepler’s problem 
whose dynamics is well known. Notice that the con¬ 
dition C > 0 is satisfied whenever the charges have 
different signs, and if the charges have equal signs it 
is required that kqiq 2 < Gmim 2 . The case C = 0 
corresponds to two free particles. The case C < 0 is 
rarely discussed in the literature, therefore we give a 
brief description of it, following the same steps as in the 
Kepler’s problem ( Goldsteinlll980l) . Notice that the sys¬ 
tem associated to C < 0 has the Kepler’s constants of 
motion, namely energy, angular momentum and those 
related to the center of mass of the system. Therefore, 
in a generic way, the movement of the two bodies oc¬ 
curs in a fixed plane, which we assume from here on. 
As first step, we make a symplectic transformation of 
coordinates, that change the vectors rq, r 2 , pi, p 2 to 
the ones 


P = 


r = r 2 - n, 

TO 1 


m i + m 2 


P 2 - 


R = 

m 2 


mi + m 2 


miri + m 2 r 2 
mi + m 2 

Pi, P = pi+p 2 . 


The upper vectors correspond to the center of mass 
of the system, whereas the lowercase ones describe the 
motion relative to the center of mass of the system. 
From the previous change of variables it is obtained 


with ro = t^-. The Hill’s region is outlined in Figure^ 

n rel 

EDITOR: PLACE FIGURE □HERE. 



Fig. 1 Hill’s region (hatched part) of the two-charged 
problem, with h* ei > 0, C < 0. 

Now we introduce polar coordinates by means 
of using the generating function S 2 {p X iPy, P, O') = 

p x pcos(9) +p y psm(0): 

x = ^ = pcos(9), 


ri = — 


m 2 


Pi = -P 


mi + m 2 
mi 


r + R, r 2 = 


mi 


m\ + m 2 


-r + R, 


mi + m 2 


P, P 2 =P 


m 2 


mi + m 2 


-P. 


By replacing the new variables and — G = |G| in the 
Hamiltonian 0 we obtain 


R*(r,R,p,P) 



IIPII 2 , \c\ 

2 M + ||r|| ’ 


( 2 ) 


where p = is the reduced mass, and M = mi + 

m 2 is the total mass of the system. The Hamiltonian 
© is separable. For the motion relative to the center 
of mass we have 


P) 


IIpII 2 

2p 


M 

INI' 


(3) 


It is not difficult to identify the Hill’s region for ©. 
This consists of all the points r on the configuration 
space such that 



where h* el > 0 is the constant of motion associated 
to ©■ Thus, the points that define this region must 
satisfy 


IMI > r 0 , 


V = =P s M0), 

Pr = ^ = Px cos (9) + p y sin(0), 

Pe = ^ = -p x psin(9) +p y pcos(9). 

With this transformation we obtain the new Hamilto¬ 
nian 

K(p,9,p r ,pe) = (pl + - f)+—, (4) 

2p \ p 2 J p 


associated to the constant of movement k* = h* el > 0 . 
In this Hamiltonian does not appear 9 , reason why the 
canonical momentum pg = l is constant, which corre¬ 
sponds to the conservation of the angular momentum. 
Replacing pg = l in ©, and considering the reduced 
energy k*, it is obtained 


k* 


Pi l 2 \C\ 

2/i 2 pp 2 p ’ 


(5) 


where is defined as the effective poten¬ 

tial V e ff(p). The constant of movement © is use¬ 
ful to obtain the phase portrait; from this we get 
p r = ±sj2p(k* — V e ff(p)). When p takes the value 

P* = §d + ^f) the canonical momentum p r is 
zero. Then, as k* grows tends to zero, and if p —> oo 
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then p r —» ±^/2 pk*. The phase portrait is shown in 
Figure [2] 

EDITOR: PLACE FIGURE [2] HERE. 




Fig. 3 Behavior of the solutions of the two-charged prob¬ 
lem with k* > 0, C < 0. 


3 The restricted charged three-body problem 


Fig. 2 Phase portrait of the two-charged problem, with 
k* > 0,C < 0. 


The ordinary differential equations obtained from 
the Hamiltonian d4]) can be integrated, as in the Ke¬ 
pler’s problem. Defining u = p , and using both con¬ 
stants of motion k* and l. we have 


thus 


d9 = - 



du 

A'lCp 2 

l2 U 


9 — 9' = — arccos 



where 9' is a constant of integration. Therefore, p and 
9 are related through the expression 


1 

P 


c [—1 + ecos (0 — 0 ')], 


( 6 ) 


where 


c = 


p\C\ 

p 


e = 



2 Pk* 
pC 2 ' 


The equation © defines a hyperbola of eccentricity e, 
whose focus is located at the origin. The hyperbola is 
sketched in Figure [3] by simplicity we have chosen 9' = 
0. Since e > 1 we have that 1 > \ = cos(9 e ) > 0 for 
some 9 e £ (0, j), thus the asymptotes of the hyperbola 
are associated to the angles 9 e , 2n — 9 e . 

EDITOR: PLACE FIGURE [3] HERE. 


In this Section we introduce the restricted charged 
three-body problem, as a limiting case of the general 
charged three-body problem. This restricted problem 
will appear when the mass of the third body tends to 
zero. 

We remember that the Hamiltonian corresponding 
to the system of bodies i = 1 , 2 ,3 of charge qi and mass 
rrii , in an inertial frame of reference, that interact under 
the gravitational and Coulomb forces on the plane, is 

U(ri,r 2 ,r 3 ,pi,p2,P3) = 


E 



E 

l<i<j<3 


where 


A ij = Gmiirij — kqiqj. 

As before, k is the Coulomb constant, G is the universal 
gravitational constant, and r* £ M 2 , p, £ R 2 , i = 1,2,3 
are the positions and momentum of the bodies 1 , 2 and 
3, respectively. 

We introduce the real parameters cq and Cij by 
means of 


qi — Oj7rq, z — 1,2, 3, 

Cij — G k(y.i OLj , z , j — 1,2,3, 

with the restriction 0:3 ^ 0. Notice that this last con¬ 
dition guarantees the existence of charge of the third 
particle, therefore a different problem from the classic 
one. 
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The differential equations associated to the problem 
of three charged bodies are given by 


Ai 2 (i~i - r 2 ) 
mi||ri - r 2 1| 3 
A 2 i(r 2 - ri) 
m 2 ||r 2 - ri || 3 
A 13 (r 3 - ri) 
m 3 ||r 3 - ri || 3 


Ai 3 (ri - r 3 ) 
mi||ri - r 3 || 3 ’ 
A 23 (r 2 - r 3 ) 
m 2 ||r 2 - r 3 || 3 ’ 
A 23 (r 3 - r 2 ) 
m 3 ||r 3 - r 2 1| 3 ’ 


(7) 


Taking m 3 —> 0 we have that q 3 —► 0, then Ai 3 —»• 0 and 
A 23 —> 0 , therefore, by 0, the position vectors ri and 
r 2 obey to one problem of two charged bodies, that is 


Ci 2 W 2 (ri - r 2 ) 
||ri — r 2 || 3 

C l i 2 m 1 (r 2 - ri) 

II r i ~ r 2 |l 3 


On the other hand, for the third position vector we have 


Ci 3 mi(r 3 - ri) C 23 m 2 (r 3 - r 2 ) 
Il r 3 — r i || 3 ||r 3 — r 2 1| 3 


The dynamics of the third particle depends on the mo¬ 
tion previously chosen for the bodies 1 and 2 , and on 
the parameters 


Ci 3 = G — kaia 3 , C 23 = G — fca 2 a 3 . 

As we saw in the previous Section, in order to have 
circular orbits for the primaries, the parameters must 
satisfy 

C7i 2 = G — ko/.\ct 2 > 0. ( 8 ) 

Remark 3.1 Without the restriction 0) on the param¬ 
eters ai and a 2 , or mi, m 2 , qi and q 2 , the primaries 
cannot follow circular trajectories, thus the circular re¬ 
stricted charged three-body problem will not be well de¬ 
fined. 

We take a circular solution for the bodies 1 and 2 (in¬ 
ertial frame), where both particles move around their 
center of mass with constant angular frequency to. We 
set ri(t) = Ci (— cos cot, — sinwf) as solution for the first 
body, and r 2 (t) = c 2 (cos ut, sinwt) for the second one; 
ci and c 2 are positive constants that satisfy the con¬ 
dition toiCi = to 2 c 2 , which corresponds to define the 
origin of the system at the center of mass of the bodies 
1 and 2. Since the bodies 1 and 2 have a circular trajec¬ 
tory, it is appropriate to take a rotating system whose 
frequency is ui. The position vectors in the inertial and 
rotating frames are related by means of the rotation 


m 


cos ojt — sin cut 


as follows: 


ri(t) = R{t)si, r 2 {t) = R(t)s 2 , r 3 = R(t)r. 


It is not difficult to see that Si and s 2 lie on the hor¬ 
izontal axis (rotating frame), that is Si = ci(— 1 , 0 ), 
s 2 = c 2 (l, 0). On the other hand, the equation that de¬ 
scribes the movement of the third body in the synodic 
frame is given by 


r = cc 2 r-|- 2 a;Jr — 


Ci 3 mi(r - si) 

||r- s i || 3 


C 23 m 2 (r- s 2 ) 
l|r-s 2 || 3 >U 


with 



1 

0 


The second order differential equation 0 can be writ¬ 
ten in terms of a pair of first order vector differential 
equations in R 2 , that is 


r = 


v, 


v = 2 co Jv + 


on 

ai 7 ’ 


where the function Q is given by 


fi = -w 2 r 2 


Ci 3 roi C 23 m 2 


Pi 


P 2 


( 10 ) 


with pi = ||r - si|| and p 2 = ||r - s 2 ||. 

It is convenient to introduce the parameters 


pi = Gmi, p 2 = Gm 2l aj 
together with 




3 = 1,2,3, 


6 = sgn(a 3 ), Pj = 1 - aja 3 , j = 1 , 2 . 

Notice that four independent physical magnitudes ex¬ 
ist: time, distance, mass and charge. We choose u> as 
unitary frequency (time), the separation between 1 and 
2 as unitary distance. Also units of charge and mass 
such that 

|CK 3 | =1) Pi + p 2 = 1- 

According to the selection of units of mass, and by the 
symmetry of the problem, it is enough to consider only 
one mass parameter belonging to (0,1/2]. For this we 
use p , where 0 < p < 1 / 2 , p 2 = p and pi = 1 — p. 
With this we have 


sinuA cos uit 


Ci 3 TOi = (1 - p)( 1 - oq<5) = (1 - p)P 1 , 
C 23 to 2 = p(l - 0 ^ 6 ) = p/3 2 - 
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Using the condition (1 — /i)ci = /xc 2 (center of mass), 
the relation Ci + c 2 = 1 (units of distance), and Gm\ = 
1 — /x, Gm 2 = /-i (units of mass), it is obtained that 
Si = (—/x, 0) and s 2 = (1 — /x, 0). The positions of 
the three bodies, in the synodic frame, are shown in 
Figure [4] 

EDITOR: PLACE FIGURE 0 HERE. 



q, 


. ^ (x.y) 


/ 1 




/ 1 

/ 

1 

/ 

1 

/ 

1 

/ 


V 

iq 2 

(T,0) 

(1-p,0) 


Fig. 4 The restricted circular charged three-body problem 
in a rotating frame with convenient units. 


restricted three-body problem of Celestial Mechanics 
(which is given by (fa, fa) = (1,1),) is in the values 
that can assume fa and fa¬ 
in the following, it is described the relation between 
the forces that the particle i = 1 , 2 , exerts on the third 
particle, according to the values that fa can assume: 

• fa < 0. The Coulomb force is repulsive and predom¬ 
inates over the gravitational one. 

• fa = 0. The gravitational and Coulomb forces have 
equal magnitude, being the second one repulsive. 

• 0 < fa < 1. The Coulomb force is repulsive and the 
gravitational force predominates. 

• fa = I. The Coulomb force is null. 

• 1 < fa. The Coulomb force is attractive. 

Denoting the coordinates of r by (x,y), and taking 
oo = 1 , we have that © is equivalent to 


x - 2 y - x = 14 , 
y + 2x - y = V y , 


( 12 ) 


The parameters that appear in (1101) have been writ¬ 
ten in terms of the new parameters fa, fa, but still 
without the restriction associated to C'i 2 . For this units, 
the condition Ci 2 > 0 is equivalent to cfacfa < 1 , which 
implies 

(fa - 1)082 - 1 ) < 1 . ( 11 ) 

The region defined by (fill) is shown in Figure [5] 
EDITOR: PLACE FIGURE [5] HERE. 


where 

/3i(l-M) fan 

P2 


V = 

with 


Pi 


(13) 


pi = V( x + P ') 2 + y 2 , p-2 = V( x - 1 + p -) 2 + y 2 - 

Considering the coordinates 

Px = x~y, p v =y + x, 



Fig. 5 The hatched part, delimited by fa = + 1, 

defines the allowed values of the parameters fa and fa. 


Remark 3.2 The difference between the circular re¬ 
stricted charged three-body problem and the circular 


we obtain that m is equivalent to the system 


x — y T Px — , p x — Ur T Py — H-x , 

y = ~X+Py = Hp y , Py = Vy — Px — —Hy, 


(14) 


whose Hamiltonian function is 

H(x,y,p x ,Py) = ]^(p 2 x+P 2 y) + (yPx ~ xp y ) - V. (15) 

The Hamiltonian function m, the potential func¬ 
tion m- and the restriction (HD. define the restricted 
charged three-body problem. 


Remark 3.3 Dionvsiou and Vaiovouloi 1 198 r \ ) consid¬ 
ered the circular restricted three-charged-body problem 
with a different parameteriza tion from the one used i n 
this work. On the other ha n d. \Radzievskn 1 195C . 1958 ). 
Kunitsvn and Tureshbaei i '198&) . \Lukvanoi 1 198A ). 


Simmom 1 198A) , studied the restricted photogravita- 


tional problem making use of a parameterization simi¬ 
lar to the one we have used in our work. In the pho- 
togravitational case it is considered, in addition to the 
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gravitational force, a repulsive (or null) force of magni¬ 
tude inversely proportional to the square of the distance, 
equivalent to the restriction pi £ (— 00 , 1 ]. 

A relevant consequence of m is that the equilib¬ 
rium solutions are confined to stay in a certain region 
of the configuration space, as we shall see in the follow¬ 
ing Section. 


4 Equilibrium solutions 

According to (TT-H) . the equilibrium solutions of the cir¬ 
cular restricted charged three-body problem must sat¬ 
isfy 

V x + x = 0, V y + y = 0, 

which implies 

x f ~ M 1 “ /*) (fj ~ f|) = °> 
yf = 0, 

where 

f = 1 _ _ P2P 

Pi P\ ' 

We can separate the solutions in two types: y ^ 0 (non 
collinear or triangular) and y = 0 (collinear). Consid¬ 
ering the first case, y ^ 0 , it is required that / = 0 , 
which implies that 



Fig. 6 Allowed region, on the parameters space (<5i, <52) = 
(di^ 3 , d^ 3 ), for the existence of triangular equilibrium solu¬ 
tions. The region is delimited by the functions 82 = <5i + 1, 

62 = <5i — 1, = 1 — Si, 5l = + 1 . 

5) - 1 



Fig. 7 Allowed region, on the configuration space, for the 
existence of triangular equilibrium solutions. 


! /3i(l — A*) P 2 P _ n 

1 o o U 

pi pi 


(16) 


and 


Pi 

pi 



(17) 


hold. The equation (fT71) requires that and have 
the same sign. If this happens, the sign must be pos¬ 
itive, so that m has solution. The corresponding so- 
lution is p 1 = pf , P 2 = P 2 ■ Therefore, it is required 
that the /3-parameters be positive, and satisfy the trian¬ 
gular inequalities pi + 1 > P 2 , P 2 + 1 > pi, Pi + P 2 > 1 . 
Besides this, the restriction CD must be fulfilled. In 
order to describe the triangular solutions on the pa¬ 
rameters space, we introduce (<3i, J 2 ) = (/3-J^ 3 , P^ 3 )', we 
will use them interchangeably. In Figures [ 6 ] and [7] it is 
shown the allowed region for the existence of triangular 
equilibrium solutions, in the parameters and configura¬ 
tion space, respectively. 

EDITOR: PLACE FIGURE © HERE. 

EDITOR: PLACE FIGURE 0 HERE. 


There are two equilibrium solutions of this type, one 
with y > 0, and another with y < 0. It is not difficult to 
see that the coordinates of these solutions are {x, y) = 
(xl,±Vl), where 

Xl = ~P + \(pI /3 - PT + !)> 

Vl = i^2(/3 ^ /3 + /3 2 2/3 ) - (P? 3 - /3 2 2/3 ) 2 - 1. 
Thus, we have proved that: 

Theorem 4.1 There are two triangular equilibrium 
solutions in the circular restricted charged three-body 
problem, whenever /3i and /3 2 belong to the region shown 
in Figured They are given by 

L4 = {x L ,y L ), L 5 = (x L ,-y L ). 


Remark 4.1 T his result was proved in th e corre ¬ 
sponding models by Kunitsvn and Turesh baeil J ’1985A) . 
Lukvanoi \19 hA . 1986) and SimmonA 1 198! i f 





























Proposition 4.1 The parameters j3i and /3 2 set the 

location of the triangular equilibrium solutions, in the 

following way: 

• p\f 3 > \J 1 + /3^ 3 . The equilibrium solution is lo¬ 
cated to the left of the body 1 (see Figure\ 8 f). 

• (3 = \J 1 + ■ The equilibrium solution is above 

or below the body 1 . 

• (3 ^ 3 < \Jl + fi ^ 3 or f3 ^ 3 < y/l + P ^ 3 ■ The equi¬ 
librium solution is located between the bodies 1 and 2 
(see Figure 0). 

• p \^ 3 = \J 1 4 - ■ The equilibrium solution is above 

or below the body 2. 

• p \^ 3 > \J 1 + p^ 3 . The equilibrium solution is lo¬ 
cated to the right of the body 2 (see Fiaure TlW . 

EDITOR: PLACE FIGURE M HERE. 

EDITOR: PLACE FIGURE M HERE. 

EDITOR: PLACE FIGURE [TO] HERE. 



Fig. 8 Equilibrium solution situated to the left of the body 
1 . 


On the other hand, for the collinear case, i.e. y = 0, 
the relative distances assume the form 

Pi = \x + p\, p 2 = \x + p-l\, 

and the roots of the equation 



Fig. 9 Equilibrium solution situated between the bodies 1 
and 2. 



Fig. 10 Equilibrium solution situated to the right of the 
body 2. 


define equilibrium solutions (by simplicity, we also use 
F(x) to denote F(x,(3 i,/ 3 2 )). The problem of the exis¬ 
tence of the roots of F depends on and P 2 , and has 
been solved in parametric form bv iLukvanov (1984), 
without the condition {pi — \){P 2 — 1) < 1. In the fol¬ 
lowing we show necessary and sufficient conditions for 
the existence of the roots of F under such restriction. 
As first step, notice that (fl8l) can be rewritten as 


F(x) = x — 


Pi (1 - p) (x + p) p 2 p(x + p- 1 ) 


Pi 


P% 


F{x,p i,/3 2 ) 


xf - /i(l - p) 




(18) 


In order to facilitate the study of the function F, it is 
convenient to separate the domain in three intervals: 


h = (-oo,-/x), I 2 = {—p, 1 — p), I 3 = {1- P,oo). 
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We have the following reductions: 

• x G I\- The relations x+p = —pi, x+p—1 = —p 2 are 

fulfilled, therefore F(x) is reduced to x + p. 

p% 

• x G h- In this case x + p = pi, x + p — 1 = —p 2 , 

then F(x) becomes x — p~ 

• x G I 3 • It is verified that x + p = pi, x + p — 1 = P 2 , 

therefore F(x) assumes the form x — 

v ' pf P2 

Note that in the three intervals it is fulfilled 

F' Or) = 1 + 20i{l - p)\ + 20 2 p\, 

Pi P 2 

where the prime denotes differentiation respect to x. 
Next, we divide the plane {0i,02) into four quadrants 
and two lines, that is 

Ri = {{0i,02) g R 2 | 0i > 0, P 2 > 0}, 

R-2 = {{Pltfa) G R 2 I 01 < 0, 02 > 0}, 

R 3 = {{01,02) G R 2 I 01 < 0, 02 < 0}, 

Ra = {{01,02) G R 2 I 01 > 0, 02 < 0}, 

R 5 = {{01,02) GR 2 I 01=0}, 

Rq = {{01,02) G R 2 | 02 = 0}, 

and define the S-regions by means of 

Sip = {( 0i,02) G R 2 | 0 < 0i < 1, 02 > 0}, 

Si ,2 = {{0i,02) G R 2 | 0i > 1, 0 < 02 < ^}, 

S 2 = {{01,02) G R 2 I 01 < 0, 02 > 

Sa,i = {{0i,02) G R 2 | 0 < 0i < 1, ^ < 02 < 0}, 

Si ,2 = {{0i, 02) G R 2 | 0i>l, 02 < 0}, 

S 5 = {{0i,02) G R 2 | 0i=O, 02 > 0}, 

5 , 6 = {(/3i,/3 2 )GR 2 I 0i >0, 02=0}. 

With this, we introduce the sets R[, i = 1, • • • , 6 , which 
correspond respectively to Ri, i = 1 , • ■ ■ , 6 , restricted 
to {0i — 1){02 — 1) < 1. Thus, we have 

R'i = Sip U Sip, 

R '2 = S 2 , 

R'a = 0 , 

R'i = Sip U Sip, 

R 5 = S 5 , 

R'e = S 6 . 

These sets define the allowed regions for {0i,02)- De¬ 
pending on the region (with exception of R' 3 ), and the 
interval, could exist one or two collinear equilibrium 
solutions. 

In the following we introduce a Theorem that re¬ 
lates different regions of the /3-parameters for which the 
collinear equilibrium solutions exist. It shall be useful 
for demonstration of subsequent results. 


Theorem 4.2 Suppose that there exist collinear 
equilibrium solutions for {0i,0^) defined by the inequal¬ 
ities 

fi{x*,p) < 0i< h{x*,p), 

gi{x*,p) < 02 < g 2 {x*,p), (19) 

hi{p) <x*< h 2 {p), 

that is, for each pair { 0 i, 02 ) there exist x such that 
F{x, 0i, 02) = 0. Then, there exist collinear equilibrium 
solutions for those { 0 i, 02 ) which satisfy 

gi{~x*, 1 - p) < 0i < g 2 {-x », 1 - p), 

fi{~x*, 1 - p) < 02 < / 2 (—z*,l- p), ( 20 ) 

—Il2(l — p) < x* < —hl{l — p). 

Proof: 

We want to show that there exist some x for the new 
/3-parameters such that this triad defines a root of F. 
Consider the change of variables x = —x, 2 * = — 2 *, 
p = 1 — p, 0i = 02, 02 = 0i- With this, the inequalities 
in (1191) take the form 

fi{~x*, 1 -p)<0 2 < / 2 (-£», 1 - p), 

gi{-x t ,,l-p)<0i<g 2 {-x t ,,l-p), ( 21 ) 

hi{\ — p)< —x » < /i 2(1 — p)- 

In the following, we also use the tilde notation for the 
variables related with x, x*, p, 0i, 02■ Notice that the 
relative distances become p\ = P 2 , P 2 = Pip the£efore 
F(x, 0i, 02) = 0 can be written as F{x,—0i,—02) = 
0. The later equation and m describe another set of 
/3-paranreters for which collinear equilibrium solutions 
exist. Removing the tilde in these relations we obtain 
(1201). 

Now we give two Theorems about the necessary and 
sufficient conditions for the existence of the collinear 
equilibrium solutions, in terms of the /3-paranreters and 
Ii, i = 1,2,3. The first Theorem deals with one sim¬ 
ple root of F, whereas the second one is related to two 
roots. The regions of existence of the collinear equilib¬ 
rium solutions are shown in Figures [TT] 1121 and fl3l 

Theorem 4.3 

1. Region R[. There exists exactly one collinear equi¬ 
librium solution for x G Ii, I 2 , -I 3 . 

2. Region R' 2 . There exists exactly one collinear equi¬ 
librium solution for x £ I 3 . 

3. Region R 4 . There exists exactly one collinear equi¬ 
librium solution for x G Ii . 
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4. Region R' 5 . There exists exactly one colhnear equi¬ 
librium solution for x £ I 3 . On the other hand, if 
p 2 > 1 holds, then there exists exactly one collinear 
equilibrium solution for x £ I\. In a similar way, if 
0 < P 2 <1 is satisfied, then there exists exactly one 
collinear equilibrium solution for x £ I 2 . 

5. Region R' 6 . There exists exactly one collinear equi¬ 
librium solution for x £ I\ . On the other hand, if 
P 1 > 1 holds, then there exists exactly one collinear 
equilibrium solution for x £ I 3 . In a similar way, if 
0 < /3i < 1 is satisfied, then there exists exactly one 
collinear equilibrium solution for i£J 2 . 

Proof: 

According to Theorem 4.2, items 2 and 4 imply 3 and 5 
respectively, therefore it is enough to demonstrate items 
1, 2 and 4. 

First item, interval I\. We notice that if a; —>• —00 
then F{x) —> —00, and if x —> —p then F{x) —> 00, 
because /3i > 0. It is clear that the function is strictly 
increasing since F'(x) > 0. The conclusion follows by 
the continuity of the function. A similar argument can 
be used to demonstrate the same item, intervals I 2 and 
I 3 , and fourth item, interval I 3 . 

Fourth item, interval I\. In this case Fix) = x+^f, 
and the function takes its maximum value at x = — p. 
Also notice that if x —> —00 then F(x) —> —00, and 
that F(x ) is strictly increasing. Therefore, the function 
could have only one root. For the existence of the root 
it is required F(—p) > 0, which implies /3 2 > 1. A 
similar argument can be used to demonstrate the same 
item, interval I 2 . 

Second item. We write the function F(x) in terms 
of p 2 = x + p — 1, as a quotient of polynomials defined 
for P 2 £ (0, 00). The equilibrium solutions are defined 
by the positive roots of the fifth-degree polynomial in 
the numerator: 

P 2 + P 2 (3 — p) + P 2 ( 3 — 2/z)+ 

+P 2 ((! - M)(l - A) - /3 2 /T)+ 


-2fi 2 pp 2 - P 2 P = 0. 

In this polynomial, the coefficients of the quintic, quar- 
tic and cubic terms are positives, whereas the sign of the 
coefficient of the quadratic term depends on the values 
of pi and p 2 , and the remain coefficients are negative. 
This implies that there is only one variation in the sign 
of the sequence of the coefficients, independently of the 
coefficient of the quadratic term, therefore we have only 
one positive root. 

In order to state the Theorem concerning two roots 
of F (possibly one root of multiplicity 2), we introduce 
a new parameter x *; it will be used to characterize the 


frontier of the region of existence of the colline ar s olu- 
i ions, on the parameters space (we proceed as Lukyanov 
( 1984 11. Once done that, we introduce 


Pl(x*,p) 


P 2 ( x *,p) 


(3a;* + p — l)(x* + p) 3 

2(1 - p) 

(3a;* + /a)(a;* + p - l) 3 
2/a 


( 22 ) 


G(x»,/a) 


P*{x*,p)P? i {x*,p) - p{{ x*,p) - P 2 {x*,p), 


and x r i(p), x r2 (p), to be defined in an implicit way. 
For fixed p, the functions x r i(p), x r2 (p) are spe¬ 
cific roots of the eighth-degree polynomial G(x»,/i) in 
a;*, for instance G{x r i{p), p) = 0. The second root 
is defined, in terms of the first one, by the expres¬ 
sion x r2 {p) = — avi(l — p). These functions satisfy 
—p < x r i(p) < — §, 5(1 ~ p) < x r2 {p) < 1 — p. At the 
Appendix we give an approximation of the roots x r \ ( p) 
and x r2 (p), using regular perturbation theory. 


The functions x r \(p) as the value of x* where the 
curve Pi = Pl(x*), p 2 = /3|(x*), —p < x* < — £ meets 
P1P2 — Pi — P2 = 0, that is G{x r i{p), p) = 0, we have 
x* £ (— p, x r i(p)) for the upper bound. 

and x r i(p), x r 2 (p) as specific roots of the eighth- 
degree polynomial G(x*, p) in x*. At the Appendix we 
give an approximation of the roots x r i(p) and x r 2 (p), 
using regular perturbation theory. 

Theorem 4.4 

1. Region R' 2 . There are at most two collinear equilib¬ 

rium solutions in each one of the intervals Ii, I 2 . A 
necessary and sufficient condition for the existence 
of the equilibrium collinear solutions with x £ Ii is 
—P*{x*) < pi < 0, P 2 (x t ) < p 2 , x* < —p. On 
the other hand, the necessary and sufficient condi¬ 
tion for x £ I 2 is P*(xf) < pi < 0, <; 

Pi( x*)-l 

p2 < P 2 ( x *), —p < x * < Xrl{p)- 

2. Region R\. There are at most two collinear equilib¬ 

rium solutions in each one of the intervals I 3 , J 2 . A 
necessary and sufficient condition for the existence 
of the equilibrium collinear solutions with x £ I 3 is 
—P 2 (x*) < p 2 < 0, P\{xf) < Pi, x* > 1 — p. On 
the other hand, the necessary and sufficient condi¬ 
tion for x £ I 2 is P 2 {x*) < P 2 < 0, *^ 2 ^*) < 

Pi < P*(x*), x r 2 (p) < X* < 1 — p. 

Remark 4.2 There are two types of collinear equilib¬ 
rium solutions involved in Theorem 4-4- Ai one case the 
root of F is simple, whereas in the other one the root 
has multiplicity 2. The roots of multiplicity 2 only hap¬ 
pen when the equalities for both Pi and p 2 are fulfilled. 











11 


For instance, 0i = —0{(xfi), 02 = /3|(x*) ; x * < ~P> * n 
the first case of R ' 2 ■ 

Proof: 

First item, interval I\. We write the function F(x) in 
terms of pi = — x — p, as a quotient of polynomials, 
defined for p\ £ (0, 00 ). The equilibrium solutions are 
defined by the positive roots of the polynomial of fifth 
degree in the numerator: 

~Pi — (2 + p)p\ — (1 + 2 p)pl+ 

+(A(i — p) + (02 — 1 )p)pI+ 

+2/3i(l — p)pi + /3i(l — p) — 0. 

Since 0 i < 0, possibly with exception of the coefficient 
of the quadratic term, the coefficients are negative, so 
it is possible to have two or none sign variations in 
the sequence of the coefficients, therefore we have three 
options: two simple roots, one root of multiplicity 2, or 
none root. This implies that F(x) is a concave function, 
so there exists a local extrema x*, that is F'(x*) = 0. 
Notice that the condition F(x*) > 0 guarantees the 
existence of at least one root. Thus, for a given p, we 
want to determine those values of 0i and 02 such that 
both F'(xfi) = 0 and F(xfi) > 0 are fulfilled. With 
this aim, as first step we will introduce a change of 
variables 0 i, 02 —> x*, fc such that x* be a root of F'{x) 
(in this case x* defines a local minimum). Next, we will 
determine the allowed values for fc; to do that we use 
F(x*) > 0 and (0i,02) G F.2- Finally, we will take into 
account the condition (0i — 1)(02 — 1) < 1 which gives 
rise to (0i,02) G R ' 2 ■ Let x* < — p the root of F'(x). 
We introduce the change of variables 

02 = ■£-{ l-x*-/r) 3 , 

2/i 

which satisfies F'(x*) = 0. The condition F(x*) > 0 
implies k > — (3x* + p), therefore 0i and 0 2 must be 
greater or equal to certain minimum values. Notice 
that k = — (3x* + p) defines the lower bound of 0i and 
0 2 , that is — /3*(x*) and /3|(x*), respectively. Besides 
this, due to ( 0 i, 02 ) G R 2 , we require 0i < 0, so at this 
point the inequalities — 0 *{x*) < 0 i < 0, 0 2 (x*) < 02 , 
x* < —p are satisfied. Since x* < —p, we have 02 > 1; 
using this and 0 \ < 0 we conclude that the inequality 
( 0 i — 1)(02 — 1) < 1 holds, therefore ( 0 i, 02 ) G R 2 - 
Second item, interval I 3 . It is consequence of the 
first item, interval I±, and Theorem 4.2 

First item, interval I 2 . We proceed as for the demon¬ 
stration of the first item, interval I\. We write the 


function F(x) in terms of p = — 1, as a quotient of 

polynomials, defined for p £ (0,oo). We focus on the 
positive roots of fifth-degree polynomial in the numer¬ 
ator: 

-0i(l - p)p 5 - 3/3i(l - p)p 4 + ((02 - 1 )p ~ 3/3i(l - 
p))p i + ((1 - 0i)(l - p) + 3/3 2P)p 2 + 3/3 2 PP + 02P = 0. 

By hypothesis 0\ < 0 and 02 > 0, therefore the 
coefficients present positive signs, possibly with ex¬ 
ception of the coefficients of the quadratic and cubic 
terms. It is possible to have two or none sign variations 
in the sequence of the coefficients, therefore we have 
three options: two simple roots, one root of multiplic¬ 
ity 2, or none root, so F(x) is a concave function. Let 
x* £ (— p, 1 — p) the root of F'(x). We use the change 
of variables 

k 

02 = 2^(1 -X* - pf, 

so F'(x*) = 0. From F(xfi) < 0 we obtain k < 
— (3 x*+p). Since (0i,02) G R 2 C R 2 we require k > 0. 
Using both previous inequalities we obtain x* < —((. 
Taking into account F(x*) < 0, it is concluded that 0\ 
must be greater or equal to a minimum value, whereas 
0 2 must be lesser or equal to a maximum value, namely 
0i(x*) and 02 (x*), respectively. Due to (0i,02) G R 2 , 
we require 0 i < Q, 0 2 > 0 , therefore we have obtained 
0t( x *) < 0i < 0, 0 < 02 < 0 2 ( x *) f° r ~P < x* < — 
The final step is to consider (0i — 1)(02 — 1 ) < 1 , 
or in an equivalent way / 3 i /32 — /3i — /32 < 0. Notice 
that the /3-paranreters in which we are interested have 
the curve (/3*(x*),/3|(x*)), —p < x* < — as upper 
bound, and /3i /32 — /3i — /?2 = 0 as lower bound. Defining 
Xri(p) as the value of x* where the curve 0i = /3*(x*), 
02 = /3|(x»), -p < x* < meets /3 i/3 2 -0i-02 = 0, 
that is G(x r i(p), p) = 0, we have x* £ (—p,x r i(p)) 
for the upper bound. With the aim of obtaining a sin¬ 
gle inequality, for the region of existence, of the equi¬ 
librium solutions, we parameterize the lower bound, 
namely 0i0 2 - 0i - 02 = 0 , using x* £ (-p,x r i(p)); 
we only must consider the part of the curve that goes 
from the point (0*(x r i(p)), 0 2 (x r i(p))) to the origin. 
The corresponding parameterization is 0i = /3*(x*), 
02 = , x * 6 (~p, x ri(p))- Therefore, the re¬ 

gion of existence of the collinear equilibrium solutions 
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is defined by the inequalities 


P*{x*) < Pi < 0 , 


/3*Qe«) 

#(*.) ~ 1 


< @2 < PU X *)’ 


—p < X* < X r \(n). 


Second item, interval I 2 . It is consequence of the first 
item, interval I 2 , and Theorem 4.2. We identify x r 2 (p) 
with —x r i(l —p), therefore —a: r i(l— p) must be a root of 
G(x*,p). According to (1221) . we have that /?*(—x*, 1 — 
p) = P 2 (x*,p) and /3*(— X*, 1 — /x) = /3|(x*,/x), which 
implies G(—a;*, 1 — /x) = G(x*,/x). From this we con¬ 
clude G(—a; r i(l — /x),l — p) = G(x r i(p), p) = 0, as 
required, for consistency. 

EDITOR: PLACE FIGURE [□] HERE. 

EDITOR: PLACE FIGURE [H HERE. 

EDITOR: PLACE FIGURE HU HERE. 


P 2 



Fig. 11 Values of the /3-parameters for which there exist 
collinear equilibrium solutions, on the interval I\. Each pair 
(/3i, P 2 ) belonging to the single hatched part allows only one 
equilibrium solution. On the other hand, the points associ¬ 
ated to the double hatched part allow two different equilib¬ 
rium solutions. Moreover, there exists only one equilibrium 
solution for the /3-parameters belonging to the continuous 
border of the double hatched region; these parameters are 
associated to a root of F with multiplicity 2. 



Fig. 12 Same description of Figure [Til on the interval I 2 . 



Fig. 13 Same description of Figure fill on the interval I 3 . 

1 /3 

3. (—n + Pi ,0) is a collinear equilibrium solution on 
the interval I3, whenever p \^ 3 — p \^ 3 = 1- 

Proof: 

Due to the similarity of the proof of these items, we will 

1 /3 

prove only item 2. Since x = —p + / 3 -p we have p\ = 
P\^ 3 ■ Notice that p 2 = plj 3 holds, since fi\^+P 2 ^ = 1. 
Thus, the condition on the relative distances p\+p 2 = 1 
is true. Finally, 

F(-p + p\ /3 ) = -p + p\ /3 - p\ /3 (l - p) + p\ /3 p 
= p(-l+p 1 1 /3 + p 1 2 /3 )=0. 


There exist collinear equilibrium solutions which 
have a very simple algebraic expression. These solu¬ 
tions are the limit of the triangular ones. 

Theorem 4.5 Assuming that (pi,P 2 ) £ R[, we have 
the following: 

l/o 

1 . (—p — Pi ,0) is a collinear equilibrium solution on 
the interval I\, whenever P 2 ^ 3 — p \^ 3 = 1. 

-l/o 

2 . (— p+ Pi ,0) is a collinear equilibrium solution on 
the interval I 2 , whenever p \^ 3 + p\f 3 = 1. 


5 Stability of the equilibrium solutions 

We know that the linear stability of the equilibrium 
solutions is determined by the eigenvalues of the matrix 

A = J RessH(x*,y*,p*,p*), 
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where (x *, y*, p* , p*) is an equilibrium solution of the 
system (fl4l) (see Mever and Hall (1992)). Then, 


A = 


( 0 

1 

1 

0 \ 

-1 

0 

0 

1 

*xx 

V xy 

0 

1 

\ Vxy 

Vyy 

-1 

0 / 


(23) 


where 


^x:r. 


fa(l-p) fap 


p\ 


P2 


+ 


+3/?i(l — p) 


(x + pf 


Pi 


3fap 


(x + p - l) 2 
p% 


V X y = 3/3i(l - p) 


(x + p)y 


Pi 


3fap 


(x + p - l )y 


p\ 


*yy ~ 


/3i(l ~ p) P 2 P 


Pi 


P% 


+3/3i(i — n)^-g + 3fan^~5- 
Pl P2 

(24) 

Using the relations pi = p{^ 3 , p 2 = , which are true 

for triangular equilibrium solutions, or the collinear 
equilibrium solutions given by Theorem 4.5, (12411 is re¬ 
duced to 


V xx = -l 


3(1 — p){x - 


3 2/3 


p) 2 + 3p(x + p-l) 2 


A 3 


V X y — 


3(1 -p){x + p)y 3p(x + p - l)y 


3 2/s 


A 3 


Vyy ~ 1 


3(1 — p)y 2 , 3 \xy 


3 2/3 


+ 


£ 13 ' 


(25) 


In order to study the characteristic polynomial associ¬ 
ated to the matrix (1231) , in each equilibrium solution of 
m, it is convenient to use fl = ^r 2 + V (see equations 
m and (1131) 1. The characteristic polynomial is 


A -|- A (4 ^xx Qyy) ^xx^yy ^ X y — 


(26) 


The triangular equilibrium solutions have been located 
using the sides of the triangle: pi, p 2 and 1. N ow we u se 
the parameterization considered by Lukyanov 198611 . 
which consists of angles 7 i, i = 1,2, and the unitary 
side; 7j is the angle between the side pt and 1, as it is 
shown in Figure [TT] 

EDITOR: PLACE FIGURE [H HERE. 

By means of using (PT5l) . the characteristic polyno¬ 
mial (l26l) becomes 


A 4 + A 2 + 9p(l — p) sin 2 (7) = 0, 


(27) 



Fig. 14 Angles 71 and 72 . 

where 7 = 71 +72, 7 S [0,7r]. Notice that the character¬ 
istic equation of the collinear equilibrium solutions, de¬ 
scribed by Theorem 4.5, is obtained by setting 7 = 0, tt. 
It is straightforward to solve (1271) using the change of 
variable u = A 2 . The eigenvalues of A are given by 


where F(p, = 1— 36p(l—p) sin 2 (7). The curve F = 0 
is defined by the equation 36 p (1 — p) sin 2 (7) = 1, with 
^-<p<\,\< sin 2 (7) < 1. Since p* = \ — 
is the smallest value of the mass parameter such that 
the function F can be zero, the associated angle must 
be 7 = f. For p £ (p*,^\, always there exist angles 
7 = 7 m , 7T - 7 P , where 7 M = arcsin( —j=L =), in a such 

way that F(p, 7) = 0 (7^ and 7r — 7^ are different from 
each other). In particular, for p = ^ the corresponding 
angles are 70 and tt — 70, where 70 = arcsin(|). 

The stability of the triangular equilibrium solutions 
is determined by the value of F; in the following The¬ 
orem we state four possible cases. Although the equi¬ 
librium solutions limit of the triangular ones (collinear 
equilibrium solutions discussed in Theorem 4.5) are not 
triangular, we have included them in the fourth item of 
following Theorem. 

Theorem 5.1 The triangular equilibrium solutions, 
as well as the solutions limit of the triangular ones, of 
the circular restricted charged three-body problem, sat¬ 
isfy: 

1. For the values of {p,”f) such that F(p ,7) < 0, the 
equilibrium solutions are unstable in the Lyapunov 
sense. 

2. For the values of (p, 7) such that F(p, r y) = 0, the 
equilibrium solutions are linearly unstable. 

3. For the values {p,"f) such that 0 < F(p, 7) < 1, the 
equilibrium solutions are linearly stable. 

4. For the values (p, 7) such that F(p,j) = 1, the equi¬ 
librium solutions are linearly unstable. 
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Proof: 

In order to cover all possible cases, observe that —8 < 
F(p, 7 ) < 1 holds for p £ (0, 5 ], 7 £ [0,7r]. 

First item. In this case the eigenvalues of the matrix 

A are A = ±\J The conclusion follows 

since the real part of the eigenvalues is different from 
zero. 

Second item. We observe that the eigenvalues of the 
matrix A are A = *, each one with multiplicity 2. 

The matrix A is non diagonalizable because for each A 
its eigenvector is 

( 2A + 14 y A 2 — 1 — V xx A 2 + 1 + V xx + XV xy A 
V D ’ D ’ D ’ )' 


where D = A 3 + A(l — V xx ) + V xy . Thus, the equilibrium 
solutions are linearly unstable. 

Third item. The eigenvalues of the matrix A are 

A = ± i\J . The solution is linearly stable, 

since all eigenvalues are pure imaginary and distinct. 

Fourth item. In this case the eigenvalues of the ma¬ 
trix A are A = 0 with multiplicity 2, and A = ±i, so A 
is non diagonalizable. The conclusion follows. Notice 
that F(p, 7 ) = 1 can only happen for 7 = 0, n. 

In Figure El are indicated the regions described by 
Th eorem 5. 1. Similar results were previously obtained 
by Lukvanov (1986). 

EDITOR: PLACE FIGURE E HERE. 



Fig. 15 Description of the values of F as function of (p, 7 ). 


5.1 Stable region in the configuration space 

In order to study the stable region in the configuration 
space, we define the restricted configuration space as 
the points (ar, y ) which meet p\ ^ 0 , pn ^ 0 , pi + 1 > P 2 , 
P 2 + 1 > Pi, Pi + P 2 > 1, ( p\ ~ 1 )(p| - 1) < 1. With 


this, we avoid collisions. Moreover, all the triangular 
equilibrium solutions have physical sense, according to 
Theorem 4.1. 

We want to show the evolution of the stable region, 
in the restricted configuration space, as p increases. As 
first step, we apply the law of cosines to the triangle of 
sides 1, pi, P 2 , and angle of interest 71 — 7 . From this 
we get 

Pi + P 2 + 2 PiP 2 cos( 7 ) = 1. (28) 


By geometry, we know that 7 = constant defines two 
arcs in the restricted configuration space, for fixed p. 
One arc satisfies y > 0, whereas the other one y < 0; we 
will refer to them as upper and lower arcs, respectively. 
With the aim of describing the arcs, we write (l28l) in 
terms of ( x,y ). After some algebraic manipulations, 
(P251) becomes 




cos(7) 
2 sin(7) 


2 


1 

4 sin 2 ( 7 ) 


(29) 


The upper (lower) arc is defined by the circumference 
(l29l) with plus(minus) sign. 

According to Theorem 5.1, and Figure El if p £ 
(0 ,p*) then the stable region is characterized by 7 £ 
( 0 , 7 r), and at p = p* appears the unstable region, as¬ 
sociated to 7 = f, so the corresponding stable region 
is defined by 7 £ (0, ^) U (f , 7 r). On the other hand, 
for fixed p £ (p* , A], the stable region is conformed by 
the points which satisfy 7 £ (0, 7 ^) U ( 7 r — 7 ^, 7 r), where 
7 u = arcsin(—J ). In the following, it is described 

6 ^//x(l-/x) 

the evolution of the stable region, in the restricted con¬ 
figuration space, according to the variation of the mass 
parameter. 

• 0 < p < p*. Exception made of the points along the 
horizontal axis, all the triangular equilibrium solu¬ 
tions are stable, since 0 < F(p, 7 ) < 1. It is outlined 
in Figure [TUI 

EDITOR: PLACE FIGURE El HERE. 



Fig. 16 Stable region (hatched part), for p € (0,/x*), in 
the restricted configuration space. 

























15 


• fi = fi*. The unstable region is conformed by the 
points on the horizontal axis, and those defined by 
'!=■§■ Except for the mentioned points, all the tri¬ 
angular equilibrium solutions are stable. It is shown 
in Figure [TT] 

EDITOR: PLACE FIGURE ED HERE. 



Fig. IT Stable region (hatched part), for fi = fi* , in the 
restricted configuration space. The dashed circumference 
indicates the unstable region associated to 7 = ^. 

• fj,* < fi < 2 - The stable region, for fixed /i, is defined 
by the points which meet 7 £ (0,7^) U (tt — 7 M , 7r), 
where 7 „ = arcsin(— . 1 ). It is outlined in Figure 


, 1 rotated ±f radians in a counterclockwise 

v 'l±cos( 7 ) 4 

sense with respect to the horizontal axis (both signs 
hold due to the dependence, of the semi-axes, on the 
harmonic function). Notice that the ellipses defined by 
7 and 7 r — 7 are related through a rotation of ^ radians. 

As we did for the restricted configuration space, we 
describe the evolution of the stable region, in the re¬ 
stricted parameters space, according to the variation of 
A 

• 0 < n < n*. With the exception of the points that 
satisfy 7 = 0 , tt, all the triangular equilibrium solu¬ 
tions are stable. It is outlined in Figure fTT)l 

EDITOR: PLACE FIGURE QJ HERE. 



EDITOR: PLACE FIGURE [TH] HERE. 


Fig. 19 Stable region (hatched part), for fi £ (0,/i*), in 
the restricted parameters space. 



Fig. 18 Stable region (hatched part), for fixed fi £ (fi*, L], 
in the restricted configuration space. The dashed arcs are 
associated to 7 = 7 ^, tt — 7 ^. 


5.2 Stable region in the parameters space 


The evolution of the stable region also can be studied 
on the parameters space (<5i, A 2 )- In analogy to what 
was done for the configuration space, we define the re¬ 
stricted parameters space as the points ( 61 , 62 ) which 
satisfy 61 ^ 0 , 62 ^ 0 , <5i + 1 > 62 , 6 2 + 1 > <5i, 
<5i + 62 > 1, (<5f — !)((£> — 1) < 1. Similar results were 
previously shown by I Si mmons (119851) , without the re¬ 
striction on /3i and fo- 

The equation <$i + 6 % + 2cos(7)<5i^2 = 1 (see (l28ll ) 
defines an ellipse with center ( 0 , 0 ), and semi-axes 


• li = fi*. The unstable region is conformed by the 
points which meet 7 = 0, n, With the exception of 
these points, all the triangular equilibrium solutions 
are stable. It is outlined in Figure l20l 

EDITOR: PLACE FIGURE [20] HERE. 



Fig. 20 Stable region (hatched part), for fi = fi *, in the 
restricted parameters space. The dashed circumference in¬ 
dicates the unstable region associated to 7 = 5 
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• /j,* < fi < 5 . The stable region, for fixed /x, is consti¬ 
tuted by the points which satisfy 7 £ ( 0 , 7 ^) 11 ( 77 — 7 M ), 
where 7 „ = arcsinf — . 1 ). It is shown in Figure 

HU 

EDITOR: PLACE FIGURE HU HERE. 



Fig. 21 Stable region (hatched part), for fixed fj, £ (/x*, |], 
in the restricted configuration space. The curves 7 = 7 ^, 7 r — 
7 m are denoted with dashed ellipses. 


6 Conclusions 

The condition (/?! — 1) (/?2 — 1) < 1 must be fulfilled 
by the parameters /3i and /? 2 , for the proper establish¬ 
ment of the planar circular restricted charged three- 
body problem. As a consequence of such inequality, 
the triangular equilibrium solutions are confined to a 
certain region of the configuration space. In a similar 
way, such restriction reduces the possible values of (3 1 
and P 2 for which collinear equilibrium solutions exist. 

As happens in the classical restricted case, the stabil¬ 
ity of the triangular solutions is guaranteed with a small 
mass parameter /i. It is interesting that, although the 
collinear equilibrium solutions (limit of the triangular 
ones) are linearly unstable, near of them appear stable 
equilibrium solutions, at least in the linear sense. 

In this work, we have considered the triangular equi¬ 
librium solutions in the plane, as well as its linear sta¬ 
bility, and the existence of the collinear equilibrium so¬ 
lutions. It would be interesting to study the stability of 
the equilibrium collinear solutions which are not limit 
of those triangular ones, as well as the equilibrium so¬ 
lutions in the space. 
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A Approximations of the roots x r \(fi) and x r 2 (n) 

In order to characterize the roots x r \, x r 2 of G(a;*,/x), we use a regular perturbation approach, considering /x 
as parameter of perturbation. Instead of work with G(a;*,/x), which is not well defined for /x = 0, we deal with 
G(cc*, /x) = 4/x(l — /x)G(cc*, /x), which is a polynomial in /x. We remember that such roots satisfy —fi < x r i < —3th 
^(1 — fi) < x r 2 < 1 — ^x. Notice that in the limit (i->0we have G(cc*,/x) —» 3x*(a;* — 1) 4 (2 + a;* + 2xl + 3cc*), 
therefore x r i —> 0 and x r 2 —> 1. For avi we have a regular perturbation problem. We approximate x r \ in a power 
series (four terms), that is 

4 

X r i — ^ ( Clif-l . 
i=l 

Equating to zero the coefficients of the powers of /x in G(x r i,/x) = 0, we get ai = — 0,2 = 0, <23 = 0, 04 = — Jy, 

therefore the required approximation becomes 

1 8 4 
= “3"“ 81^ • 

One the other hand, x r 2 cannot be handled directly by a regular perturbation approach. Nevertheless, by means 
of the change of variable 2V2 = 1 + ez, where e = /x 1 / 4 , we obtain a regular problem with 2 as variable, and e as 
perturbation parameter. Using 

3 

Z = Y1 bi£l 

i=0 

in G(1 + e2, /x) = 0 we get b 0 = h = b 2 = b 3 = |||, thus we obtain the approximation 


x r 2 = 1 — 



11 .4 , 67 497 

36x/3 ^ 2 + 864^/12^ 4 486^ 


In this case, we had two possible values for b 0 , one positive, the other one negative. In order to satisfy the 
inequality x r 2 < 1 — /x, the negative vale was chosen. 
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